Accurate measurement of the D° — D° mixing parameters 
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We propose a new method to determine the mass and width differences of the two D meson 
mass-eigenstates as well as the CP violating parameters associated with D° — D° mixing. We show 
that an accurate measurement of all the mixing parameters is possible for an arbitrary CP violating 
phase, by combining observables from a time dependent study of D decays to a doubly Cabibbo 
suppressed mode with information from a CP eigenstate. As an example we consider D° — > K*°tt° 

and K s tv°. We also show that decays to the 



decays where the K*° is reconstructed in both K + n 
CP eigenstate D -> K + K~ together with D -> K + tv~ 
parameters. A combined analysis using D° -» K*°tv° and D 
the ambiguity in the determination of parameters. 



decays can be used to extract all the mixing 
K + K~ can also be used to reduce 



PACS numbers: 14.40.Lb, 11.30.Er, 13.25.Fc, 12.60.-i 



Evidence for mixing in the neutral D meson system has 
recently been reported [l], 0, [1] by the Belle and BaBar 
collaborations. These experiments find non- vanishing 
width and mass differences between the two neutral D 
mass eigenstates assuming negligible CP violation. In 
this letter we propose a method to determine all the mix- 
ing parameters accurately allowing for arbitrary CP vio- 
lation. 

Within the Standard Model CP violation in the D 
system is negligible. Hence observation of CP violation 
would be a good sig nal for New Physics (NP) 0L While 
no CP violation has been seen in D — D mixing 3[ , with 
the current precision large possible NP contributions are 
not ruled out. 

We show that using the doubly Cabibbo suppressed 
(DCS) mode D — > A'*°7r° and its conjugate modes, we 
can solve for all the D — D mixing parameters. This 
is possible if the K*°/K*° is reconstructed both in the 



self tagging K 



mode and in the CP eigenstate K s n 



mode. While the CP eigenstates D — > K + K~ cannot 
alone be used to determine all the mixing parameters, wc 
demonstrate that minimal additional information from 
DCS modes allows determination of all parameters. This 
approach may provide the optimal method to determine 
all the parameters with current data. In both these cases, 
the parameters can be determined accurately even in the 
limit of a small or vanishing CP violating mixing phase <j>. 
It has recently been proposed to use the singly Cabibbo 
suppressed (SCS) D — > K*K modes to determine the 
mixing parameters [H, @]. However, if <p is zero, these 
methods would be feasible only if the strong phase in- 
volved is measured elsewhere. The strong phase can be 
measured using a Dalitz plot analysis [7|. In the ab- 
sence of the strong phase information, these modes can- 
not be used to determine the mixing parameters accu- 
rately when is small, since they can only be expressed 
as ratio of small quantities. 



Our study of the various modes allows us to conclude 
that in the limit of small 0, an accurate measurement of 
all mixing parameters is possible only if the method also 
allows the determination of the parameters in the case 
<t> = 0. While mixing parameters can be determined us- 
ing decays to SCS non-CP eigenstates alone, an accurate 
measurement of mixing parameters in the limit of small 
4> is possible, only by adding information from decays to 
CP-eigenstates or if the strong phase is measured inde- 
pendently elsewhere [1, 0, 0, H 0, El ■ The D -» K* tt° 
modes are an example where it is possible to measure all 
the mixing parameters and the strong phase using only 
related final states, thereby reducing systematic errors. 
The methods discussed in this letter do not have sys- 
tematic errors associated with the parameterization of 
the resonant content of the Dalitz plot 0] and hence are 
model-independent. 

The neutral D mass eigenstates are related to the weak 
eigenstates by, \D 1)2 ) = p\D°) ± q\D a ) . The mass and 
width differences of these eigenstates are popularly writ- 
ten ll| in terms of the dimensionless variables, 
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where T is the average of the widths of the two mass 
eigenstates. If the magnitude of q/p differs from unity 
and/or the weak phase <f> = &rg(q/p) is nonvanishing, 
this would signal CP violation. We consider mixing to 
be the only source of CP violation and assume that the 
decay amplitudes themselves have no weak phase 12[. 

In the limit i « 1, i/ « 1 and Tt <C 1, the time 
dependent decay rates for a D° decaying to a final state 
/ and D° — > / and D° — > / have the form: 



\A(D°(t) ^/)| 2 = 
\A(D°(t) ^/)| 2 = 
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We first consider the DCS mode D° — > K*°-k° and its 
conjugate mode D a -> ^*°7r°, with the K*°/K*° recon- 
structed in the self-tagging K ± n T modes. The coefficient 
functions of the constant, linear and quadratic terms in 
(r t) in the time dependent decay rates are given by, 
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where, 



(x cos (5 A -* X + y sin 5^ ) , 
( y cos <5k» t — x sin 6^ ) , 



(3) 
(4) 
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A(D° 



K*°n°) and the ratio of the DCS 



with A K 
to CF amplitude defined as 

A(D° -> K*°ir°) _ A(D° -> if *°7r°) 



The amplitude | A^ | can easily be measured using the 
time integrated rate for the Cabibbo favored (CF) mode 
D° — > K*°ir° which is given by, 



\A(D v {t) K*\ u )\ 2 dt « \A 



(10) 



where terms of the order of x 2 or y 2 and a; or y 
are neglected compared to unity, as these are expected 
to be 0(1O~ 4 ) or less. The ratio can be determined 
using Eqs. ((3|) and (fT0|) . The observables Zk% and 
also readily determine \q/p\ and a; 2 + y 2 to be: 



l-l 4 

fy 2 



V Zk%Z k 

: i r?~ 
4*, 



(11) 

(12) 



The two linear terms in the time dependent DCS decay 
rates and may be re-expressed in terms of two 
more convenient observables Yi +) and Yl~ as follows 
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2r /rtr |^„| 2 |g| |p| 
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(13) 
(14) 



Note that the observable Y^ ) may be difficult to measure 
in the small d> limit. 



The K*°/K*° in the final state could also have been 
reconstructed in the neutral K s tt° mode, resulting in an 
additional observable. A unique feature of the final state 
K s tt°tt is that it includes contributions from both K*°tt° 
as well as K*°ir° states; the amplitude for this final state 
is thus a sum of the CF and DCS amplitudes, 



\A K 



\A(D° ^ K s tt\°)\ 2 

| Atn, | 2 ( 1 + r 2 ^ - 2 cos 4*, 



(15) 



Since the decay mode involves two neutral pions it will 
not be easy to perform a time dependent study. Hence, 
we consider only the time integrated decay rate for this 
mode. The amplitudes A(D° -> K s ir°Tr°) and A(D° -> 
K s ir°Tr°) are equal since K s tt°tt is a CP eigenstate. 
Hence, the time integrated decay rate for D° — > if s 7T 7r 
is given by: 



\A{D°(t) -> K s n\°)\ 2 dt 
'■ |A<T7r7r| 2 |"l + -(y cos<j) — x sin 0)1 

p 

! 1^4^*^ | 2 riH — (y cos0 — x sine!))— 2r x % cos 5^1 ,(16) 
p 



where, terms of order x 2 , y 2 and r K \ x, r K \ y as well as 
rz* are once again neglected compared to unity. 



Using Eqs. ([13]) and (|14|) , one obtains the following 
solutions for tan 2 tf> and y' 2 ,^ : 
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where, fa = P - XX + XX ■ The ambiguity in the 
solutions of the quadratic equations in tan 2 <p and y' 2 % is 
fixed by the correct limiting solution in the = limit. 
Further, expressing cos in terms of x, y and x'^, y'^, 
Eq. (|16[) may be rewritten as a quadratic equation in x/y, 
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allowing x/y to be solved with a four-fold ambiguity, a; 
and y can thus be individually determined using Eq. (|12| . 
The solution obtained is finite even if = 0, with a cor- 
rection term of order Y.i~ ) . Hence an accurate estimation 
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is possible even if cj> is tiny. We show below that the am- 
biguity in x/y can be reduced if information from K + K~ 
modes is added as well. 

We next consider the time dependent decay of a D me- 
son to a singly Cabibbo suppressed (SCS) CP eigenstate 
such as D — ► K + K~ or D — > tt + it~. To be specific, we 
will consider only the K + K~ final state, but the con- 
clusions can be straightforwardly applied to any other 
SCS-CP eigenstate. For this final state, the strong phase 
is identically zero; and hence, the coefficients of the con- 
stant and linear terms in (Tt), defined using the time 
dependent decay in Eq.Q]) reduce to the simple form: 
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Yx K ' . We once again examine 



where T KK = f 
in detail the solution for the case of small 4>. If cj) is 
small, the measured value of Ykk 1 will be small. The 
above solutions can then be written as a series in YSz' : 
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and therefore x 2 /y 2 is finite even for small tfi. 

As mentioned earlier, if information from K + K~ 
modes is added to that from the K*ir modes, a reduction 
in ambiguity is possible. If the observables Y^p and Y,^ 
are used, then cos 8 K \ can be obtained purely in terms 
of observables directly from Eq. p^|) . Knowing cos 5k\ , 
Eqs. (H3J), dHJ) and jM]), (J25J can be used to get, 



Unlike the DCS modes where the term quadratic in Tt 
is enhanced by the ratio of CF to DCS rates, in the SCS 
modes all time dependent terms are of the same order 
in sin# c , hence quadratic and higher terms in cannot 
be extracted. Assuming \q/p\ « 1 and = 0, the linear 
term in T t can directly measure y as has been done in 
Ref. However, the time dependent study of only the 
SCS CP eigenstates does not allow x to be determined, 
even in the limit \q/p\ ~ 1 and = 0. 

We will show that if we also include in this analysis 
the quadratic terms in (r t) from the time dependent 
decay rates of DCS modes such as Kit, all the mixing 
parameters can be solved without approximation. For 
D° — > K + tt~ and D° — > K~tt + , the coefficient func- 
tions of the quadratic terms in (r t) will be analogous to 
those for the K *tt mode given in Eq. ([6|). Hence, the cor- 
responding observables Z Klr and Z Kl , readily determine 
\q/p\ and f 2 = x 2 + y 2 . Alternatively, \q/p\ and f 2 could 
be measured using time integrated wrong sign relative 
to right sign scmilcptonic decay rates. Having obtained 
\q/p\ and / 2 , and x/y can easily be determined from 
D -> K+R-. Using Eqs. |[2"T] ) - ([25] ) . which can be re- 
expressed as, 
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solution for x 2 /y 2 and can be straightforwardly written 
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Combining this with the sum x 2 + y 2 , x 2 and y 2 can be 
individually determined. Further, Eqs. ([25)) and (fT4)l can 
be used to obtain, 
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which helps in reducing the ambiguities in x and y from 
four- fold to two-fold. 

Recently a method was proposed [B| to determine all 
the mixing parameters using the D — > K* K modes. This 
mode is singly Cabibbo suppressed (SCS) but unlike the 
K + K~ mode it is not a CP eigenstate. Hence one can 
study time dependence in all the four modes: D° — > 
K^KT and D° -> K^K*. The coefficients of the 
constant and linear terms in (r t) may be written as: 
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where, A K . K = A(D° -> K*+R-) = A(D° K*~ K+) 
and v k*k is defined as 
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It may also be noted that J/^._ K+ and x'.^_,^ are different 



U K*-K+ 



from y' I( t+ K - and x'^^ and are defined as: 



x 



/ 



(x cos <W ± y sin <W ), 

(y cos Stat T a: sin 5 K * K ). (28) 



One may conclude that the six obscrvablcs in Eqs. ((27|) 
can be used to evaluate the six parameters |Ak* k | , r 2 ^, 
x, y, 4> and 5 K * K assuming the value of \q/p\ from else- 
where. However, note that if the mixing phase = 0, 
then, the number of observables reduces to four (since 
now, IffpiW- = |i?| 2 *W- and |g| 2 iWf = \p\ 2 Y^-^) 
and a solution of all the five parameters is not possi- 
ble without some additional information. Moreover, for 
small but nonvanishing (j> the solution for the ratio x 2 jy 2 
will be inaccurate as it will depend on ratio of two very 
small observables 



To sec this, let us define, Y^_ 



sm < 



and Y£> 



y(+) 
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x' lf ^_ K + sin </>, which can all be determined in 



terms of observables using Eqs. 
that: 



(|27)1 . It then, it follows 
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It is clear that the RHS involves the ratio of two small 
quantities when <f> is small. It is easy to see that this situ- 
ation is easily alleviated if 8k> k is measured elsewhere Q ■ 
In fact, the knowledge of 5^ not only allows the ad- 
ditional determination of \q/p\ @, but also enables an 
accurate measurement of mixing parameters. 

We now estimate the values of the mixing parame- 
ters that can be obtained using the current data for 
D — > K + K~ /-k + -k~ ,1] and the world average for 
x? + y 2 [l5|]. Assuming \q/p\ = 1, we obtain Y^P = 
0.0131 ± 0.0041, = -0.0001 ± 0.0034 and f 2 = 

0.00042 ±0.00022, resulting in |x| = (1.57 ± 0.56) x 10" 2 , 
\y\ = (1.31 ± 0.41) x 10 -2 and value (up to ambiguities) 
of <fr = ±(0.36 ± 12.36)°. We emphasize that our method 
allows the determination of mixing parameters even for 
\q/p\ ^ 1; the choice \q/p\ = 1 has been made here, only 
due to lack of complete tabulated data. 

An estimate of the precision to which the mixing pa- 
rameters can be measured, using the D — » K*tt modes, 
requires the number of reconstructed D — > K*°ir° — > 
K + n~ir° events. While a branching fraction for this 
mode has not yet been reported, about 500 events (in 



K+ 



7T 7T 



havc 



230/6" x ) for the mode D° -> K 
been observed We present our estimates using two 
representative values for the ratio of DCS modes: 



B(D° 



B(D° -> K* 



K + TT TT 



= (0.4,1.2). 



These values are chosen to be of the order 0.85, the mea- 



sured 1J] ratio of corresponding CF branching fractions. 



With an integrated luminosity of 1 ab 1 at an e + e B 



factory we expect about (4000, 12000) D° -> . 
A' + 7T _ 7r° events. Interpolating the errors in D 



and assuming 5 



0, the approximate errors on 



\x\ and \y\ are expected to be (4.7,2.7) x 10 and 
(8.9,5.2) x 10~ 3 , respectively. 

We have proposed a new method to determine the 
D — D° mixing parameters x, y, \q/p\ and <f> for arbi- 
trary values of <p- The doubly Cabibbo suppressed mode 
D° — > A*°7T reconstructed in two final states (K + tt~tt° 
and A s 7r°7r°) enables the determination of all the mix- 
ing parameters. For the K s ir°ir° mode, only time inte- 
grated measurements are used, while for the A + 7r~7r° 
mode time dependent measurements are required. We 
also show that decays to the CP eigenstate D — > K + K~ 
together with D — > K + ir~ can be used to extract all 
the mixing parameters. By combining measurements of 
D — > A'*°7T with results on D — > K + K~ one can reduce 
the number of ambiguous solutions for mixing parame- 
ters. We estimate that \x\, \y\ and <j) can be measured 
with precision of order 0.6 x 10~ 2 , 0.4 x 10~ 2 and 12° re- 
spectively, using data available at present. It should be 
possible to determine |x|, \y\ to order 7 x 10~ 4 , 4 x 10~ 4 
respectively and <\> to about 1° at a Super-B factory with 
an integrated luminosity of 50 a6 -1 [16l |. 
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